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Abstract 

We calculate the contribution to inclusive high tranverse momentum (pt) charm pro- 
duction at HERA from the excitation of charm in the photon. At large values oipx the 
results of such a calculation, in the structure function language, will be more reliable 
as it sums the large logs, log(py/m^), as opposed to calculating the contribution of the 
2—5-3 subprocess in fixed order of perturbation theory. We find that this contribu- 
tion is large and comparable to the contribution from 75 fusion production of charm. 
Suitable cuts on the rapidity of the 'away-side' large pT jet allow a very neat separa- 
tion between the contributions from the excitation process and from pair-production. 
We further find that including this excitation contribution we can reproduce the mea- 
sured inclusive D* and cross-sections measured by the ZEUS and HI collaborations 
respectively, in a LO calculation. 
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Measurements of F2 in 7*7 scattering at the e"'"e~ colliders PEP, PETRA, TRISTAN 
and LEP have by now yielded a lot of information on the parton content of the photon 
over a wide range of x and Q^. However, these measurements give direct information only 
about the quark content of the photon. The gluon density g^{x,Q'^) is poorly determined 
as it affects F2 only through the QCD evolution equations. At the current values of 
the charm quark contribution to F2 is approximated by the quark-parton-model (QPM) 
matrix elements for the process 77* — > cc and 'j*g'^ —>■ cc. Through the latter process, 
the effective charm content of the photon becomes sensitive to g'^{x,Q'^). At larger values 
of Q^, c^{x, Q^) computed using the massive Altarelli-Parisi (AP) evoultion equations, is 
considerably different from the pure QPM predictions A study of the charm content 
of the photon might also help shed some light on the correct treatment of a heavy parton 
inside a target. The various different available parametrisations of q'^{x,Q'^) and g'^{x,Q'^) 

treat the charm density c'lx, Q^) with varying amount of rigour and care. It is therefore 
interesting to take a phenomenological approach and think of measurements which will probe 
c'{x,Q'^) directly and hence perhaps also yield information about g"'{x,Q'^). 

One possibility is to study production of single charm in ep collisions via the excitation 
processes (the subprocesses being c' + c + q and c' + g c + g, here we neglect the 
contribution coming from charm in the proton) shown in fig. ??. This will give rise to a 
single high— charm particle whose transverse momentum is balanced by a light q/ g jet. 
Of course the use of structure functions to compute this process is meaningful only for large 
values of the px of the charm quark. Admittedly for lower values of pt the more reliable 
computation will be that of the 2 — 3 subprocesses (some of which are shown in fig. 2), but 
at larger values of px the structure function language sums up the large \og{p\lm1) terms 
and hence is more accurate. Another contribution to the inclusive charm signal comes from 
cc pair production, via the 'direct' -^g fusion subprocess as well as the 'resolved' processes, 
where the balancing high-pj^ jet is the c (c) quark jet. 

The excitation contribution of diagrams in fig. ?? is given, in the Weizsacker- Williams 
(WW) approximation, by 
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where /^|e, /c|7 and fp^\p represent the flux factors of the 7 in the electron, charm in the 
photon and parton Pi in proton respectively; 
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where z^^f^^^^^^-^ correspond to the experimental cuts on the outgoing electron (or equivalently 
the 7 energy) and z^^^, x™" and x™*" correspond to the kinematic limits on the different 
momentum fractions. The direct contribution to cc pair production is similarly given by 
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where fg/p represents the gluon flux in the proton, and the 'resolved' contribution is given 
by an expression similar to Eq.(P where contributions from all the various subprocesses 
involving all the different partons in the photon are to be included. We use LO expressions 
for all the 2 — > 2 subprocess cross-sections 0. 

The virtuality —P^ of the exchanged photon in figs. ?? and 2 can, in principle, affect 
its parton content Q]. For the results presented here we impose the requirement < 
0.01 GeV^ and 0.25 < z < 0.70, following the cuts used in the experimental study of the 
photoproduction of jets 0. This implies that our expression for the photon flux factor is 
given by 

a ,^ , /O.Ol GeV^X a\ - z 

+ --— , (3) 

where 
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As a result of the cut on P^ we can neglect the effect of the virtuality of the 7 on its parton 
content. 

The results of our computations are shown in fig. ?? for various photon structure function 
parametrisations DG |^, LAC0], WHIT0 and different proton structure functions p|; we 
use Aqcd = 0.4 GeV for the DG and WHIT parametrizations, and 0.2 GeV for LAC. We see 
from this figure that the excitation cross-sections are indeed comparable to the cc production 
cross-sections. This implies that while the 'resolved' contribution to cc pair production is 
small (as shown by the long dashed line) for these large values oipT, the inclusive charm signal 
still has a considerable 'resolved' component due to the excitation contribution. Though we 
do not show them separately here, the contribution to a'^^'^ coming from a gluon in the initial 
state dominates over most of the pt range. Apart from the DG parametrisation for which 
excitation contributions are about a factor 2 higher than the rest, the charm excitation 
cross-section at HERA seems to be fairly independent of both the photon and the proton 
structure function parametrisations. Since the DG parametrisation has c'^ = u'^ it definitely 
overestimates the charm excitation and hence this part of the result is easily understood. 
In principle, the other parametrisations of c'*'(x, Q"^) do also look quite different, both in the 
small and large region, but the effective c— quark content of the electron, which involves 
the convolution of this with the WW function is very similar in the end; this is reflected in 
the similarity of the predictions using the LAC and WHIT parametrizations for c'''(x,Q^). 
This can be looked upon as a positive point in that the size of the excitation contribution to 
the inclusive charm signal, at large pt, can be estimated quite reliably. However, it is also 
clear that one needs to devise kinematic cuts to separate the excitation contribution to the 
inclusive charm signal from that due to cc pair production. (This is also necessary if direct 
cc pair production is to be used to study the gluon density in the proton.) 

To this end, we next study the kinematic distributions of the decay muon which is used 
to tag the charm in the final state and also that of the balancing ('away-side') jet. In this 
calculation we include fragmentation of the c— quark into a charmed meson a la Peterson 



fragmentation function with the parameter e as given in ref. |I0[, and use the value 0.1 
for the semileptonic branching ratio of the charmed hadrons. Since, for our pt cut, cc 
production is dominated by the 'direct' process, the real kinematical difference between 
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the excitation and cc contributions to the charm signal comes from the fact that in the 
excitation process only a fraction of the 7 energy is available for the subprocess, whereas 
in the 'direct' process all of it goes into the subprocess. As a result, the direct process on 
the whole receives contributions from smaller Xp values than the excitation process does. 
Hence the c jet in the cc case will have much more negative rapidity than the 'light' (q/g) 
jet in the excitation process (the proton direction is taken as positive z axis); this is very 
similar to the corresponding situation with the photoproduction of jets 0. On the other 
hand, the rapidity distributions of the charm quarks produced in both the excitation and 
pair production processes are very similar and hence those of the decay muons also. The 
kinematic distribution in pj. and therefore are very similar for both contributions. Fig. 
4 shows the rapidity distribution of the jet balancing the large pt charm, with the WHITS 
parametrisations of q^i^x^Q"^) and g^i^x^Q"^) and MRSD- for the proton structure function. 
As we can see very clearly from the figure, a cut on yj^t < 0.5 can neatly separate the 
excitation and the cc contribution from each other. The rates presented in the figure include 
the semileptonic branching ratio of the charm meson. It should also be mentioned that these 
distributions do have some sensitivity to , but only for negative values of Ujet where the 
signal is dominated by the cc contribution. The figure also tells us that the signal is healthy 
even after these cuts and hence is measurable. For a clear signal one will have to make an 
additional cut on pj. as well but that will affect both the excitation and the pair production 
contribution similarly. 

Recently both HI and ZEUS have reported measurements of inclusive charm production 



at HERA |TT1. ZEUS reports D* production with Pt{D*) > 1.7 GeV and \r]{D*)\ < 1.5 
where r] is the pseudorapidity with < A GeV^ and 0.15 < 2; < 0.86, whereas HI reports 
observation of events with a hard muon with pi^ > 1.5 GeV and 30° < 6{^) < 130°. The 
ZEUS analysis then uses this measured cross-section to estimate the 'total' cc cross-section 
by extrapolating it outside the measured region and then compare the value so obtained 
with the QCD NLO calculations. We attempted instead to reproduce the cross-sections 
measured by ZEUS and HI by using our LO QCD calculations. Since the pr cut and nic are 
comparable it is not clear whether factorisation of the production and fragmentation of the 
charm quark is such a good approximation. On the other hand if we include the fragmenta- 
tion of the c— quark in the final state then we must include the excitation contribution which 
corresponds to the fragmentation of the initial state photon into charm. We therefore run 
our Monte Carlo with two different options: In one case (A) we fragment the charm using 
the Peterson fragmentation function and include the excitation contributions whereas in the 
other case (B) we do not include the fragmentation of the final state c quark and drop the 
excitation contributions as well. Eventually detailed comparisons with measured transverse 
momentum and rapidity distributions should reveal which discription is more appropriate. 
When comparing with ZEUS results we include a factor of 0.26 which is the probability of 
a charm quark to fragment into a charged D* meson. 

Table |I] gives a summary of our calculations along with the results reported by the 
two experimental groups. We find that for the ZEUS data, in case A the results become 
less sensitive to the low- a; behaviour of the proton structure function and the excitation 
contribution actually dominates. In this case both the MRSD- and MRSDO for the partons 
in the proton and WHIT5 or LACl partons in the photon reproduce the data whereas DC 
predicts too big a cross-section. On the other hand, in case B, the results are more sensitive 
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Table 1: The D* and /i cross-sections measured at HERA by ZEUS and HI, compared 
with the LO predictions discussed in the text. 



ZEUS {D*) 


HI (/.) 


Data: 32 ± 


nb 


Data: 2.03 ± 0.43 ± 0.7 nb 


A: Era} 


5 and c' included 


A: Erag and c' included 


p str.fn. 


7 str. fn. 


a (nb) 


p str.fn. 


7 str. fn. 


cr (nb) 


MRSD- 


LACl 


27.4 


MRSD- 


LACl 


l.S 


MRSD- 


WHITS 


26.0 


MRSD- 


WHITS 


l.S 


MRSDO 


WHITS 


21.9 


MRSDO 


WHITS 


1.48 


MRSD- 


DG 


64 


MRSD- 


DG 


2.9 


B: No frag 


and c' not included 


B: No frag and not included 


p str.fn. 


7 str. fn. 


a (nb) 


p str.fn. 


7 str. fn. 


a (nb) 


MRSD- 


LACl 


26.8 


MRSD- 


LACl 


1.8 


MRSDO 


WHITS 


1S.4 


MRSDO 


WHITS 


l.S 


MRSD- 


WHITS 


37.0 


MRSD- 


WHITS 


2.4 


MRSD- 


DG 


33.4 


MRSD- 


DG 


2.3 



to the low-x behaviour of the parton densities and we find that we can reproduce the cross- 
section only for a steeply rising gluon density. In this case all the photonic parton densities 
combined with MRSD- for the proton are acceptable whereas MRSDO gives answers smaller 
by a factor 2 for all reasonable choices of the momentum scale as well as the photonic parton 
densities. Eor the HI sample, the pj. cut means that the produced charm quark has much 
higher pt than for ZEUS. Again, inclusion of fragmentation reduces the sensitivity to the 
low-x behaviour of the gluon in the proton and results for various combinations of partons 
in the photon and proton are almost the same. As before, the DG parametrisation predicts 
a large cross-section, 2.9 nb, but it is not inconsistent with the data. 

In conclusion we have studied the contribution to the inclusive charm signal from the 
excitation of charm in the photon. We find the rates to be comparable to the contribution 
coming from cc pair production. Due to the convolution with /^|e the sensitivity of a^^'^ of 
Eq.(|l]) to the region of large x^, where the various parametrisations for c'{x, Q"^) differ most, 
is reduced. This means that this 'resolved' background to the large pxi inclusive charm signal 
coming from the 'direct' process can be predicted quite reliably. Making a cut on the 'away- 
side' jet-rapidity allows a separation of these two contributions. Note that charm excitation 
events should contain a second 'spectator' charm (anti-)quark in the photon remnant, which 
has a large, negative rapidity (close to the electron beam direction). In contrast, cc pair 
production events should contain a second high— charm (anti-) quark. In both cases the 
presence of a second charmed particle should be visible in some fraction of the events; this 
could be used as a cross-check of the relative sizes of the two contributions to the inclusive 
charm signal. Even after all the cuts and folding of the cross-sections with various branching 
fractions, the rates for inclusive charm production are large and easily measurable. Eurther, a 
LO computation including the excitation contributions gives results comparable to the recent 
measurements of the inclusive charm signal by ZEUS and HI. On the theoretical side, it 
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would be interesting to compare our results with explicit 2 — > 3 subprocess calulations and 
see at what values of pr do they match. 
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